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f-^^ _ Abstract 

O , . Thermodynamical properties of nuclear matter are studied in the frame- 

work of an effective many-body field theory at finite temperature, considering 
the Sommerfeld approximation. We perform the calculations by using the 
nonlinear Boguta and Bodmer model, extended by the inclusion of the fun- 
damental baryon octet and leptonic degrees of freedom. Trapped neutrinos 
are also included in order to describe protoneutron star properties through 
the integration of the Tolman-Oppenheimer-Volkoff equations, from which we 
k>( ^ obtain, beyond the standard relations for the masses and radii of protoneu- 

j_] . tron stars as functions of the central density, new results of these quantities 

C^ ' as functions of temperature. Our predictions include: the determination of an 

absolute value for the limiting mass of protoneutron stars; new structural as- 
pects on the nuclear matter phase transition via the behaviour of the specific 
heat and, through the inclusion of quark degrees of freedom, the properties of 
a hadron-quark phase transition and hybrid protoneutron stars. 

Keywords: Quantuni-Hadrodynamics, Finite Temperature, Protoneutron Star 

1. Introduction 

In the last few decades, studies on the evolution of compact stars became a 
central issue related to theoretical and experimental research on the equation of 



*E-mail: gfm@if.ufrgs.br 



2 G.F. Marranghello, C.A.Z. Vasconcellos, M. Dillig and J. A. de Freitas Pacheco 

state (EOS) of dense matter. In 1967, the first pulsar was observed"'^ and, based 
on characteristic observational features, this object was identified as a rotating 
neutron star. Since this direct evidence for the existence of neutron stars, nuclear 
models have been widely employed in the description of the internal structure, 
composition, dynamics and evolution of these massive compact stars. As under 
the pull of gravity, the energy density in the core of compact stars is thought to 
approach or even exceed the critical value of 1 GeV/ fm? , i.e. more than 6 times 
the density of nuclear matter, their structure depends sensitively on the equation 
of state for very dense matter. As the 

EOS describes how the energy density and pressure vary with density and tem- 
perature, 

it should be also able to describe different phases, including gaseous nuclear 
matter and liquid nuclei structure up to the deconfinement transition. 

The EOS provided by model calculations and the Tolman-Oppenheimer-Volkoff 
(TOV) equations^'"^ characterize the structure of the star, whose properties like 
mass, radius, crust extent, moment of inertia are then susceptible to be compared 
to observations'*. Moreover, these models should also provide the description of 
relevant dynamical properties of compact stars as the rotational period, the emission 
of neutrinos and gravitational waves as well as more fundamental aspects involving 
the evolution, internal composition and structure of these stellar objects. 

Concerning neutron stars, nuclear matter models should be able to account for 
such objects at least as massive as the most massive observed pulsar, representing 
an important constraint on any proposed EOS for dense matter. Presently, the 
best mass determination corresponds to the binary pulsar PSR 1913+16 (M = 
I.444M0). Other estimates based on quasi-periodic oscillations (QPOs), observed in 
the X-ray emission of low mass binaries (LMXBs) suggest values up to 2.0 — 2.2M0, 
but these results are quite uncertain since they are model-dependent. More recently, 
from new radial velocity data, the mass of the neutron star associated to the X- 
ray source in the system Vela X-l"'' was recalculated with a larger accuracy than 
precedent estimates: the resulting value is 1.86 ± O.I6M0. For a given family of 
stars, stiffer EOS predict higher limiting mass and larger radii (i? > 10— 12 km). In 
particular, recent hadrodynamics models point to EOS stiffer than those expected 
on the basis of a non-relativistic approach, in better agreement with Vela X-1 results. 

In the theoretical treatment of properties of compact stars, the relativistic phe- 
nomenological approach developed by Walecka®, quantum hadrodynamics (QHD), 
represents one of the most important approaches to the highly nonlinear behav- 
ior of QCD at the hadronic energy scales. QHD model is a relativistic quantum 
field theory based on a local lagrangian density which uses the nucleon and two 
meson fields, the scalar-isoscalar, attractive a and the vector-isoscalar, repulsive w, 
as the relevant effective degrees of freedom. This model provides a thermodynam- 
ically consistent theoretical framework for the description of bulk static properties 
of strong interacting many-body nuclear systems under extreme conditions. The 
nonlinear model of Boguta and Bodmer (BB)^ was developed to improve the de- 
scription of static properties of nuclear matter, relative to the original Walecka 
model. With the introduction of cubic and quartic scalar meson self-interactions in 
the lagrangian density, the model has two additional parameters which allow and 
provide sufficient fiexibility to reproduce at saturation density, currently accepted 
values for the compression modulus of symmetric nuclear matter and the nucleon 
effective mass. 

In this work we investigate static properties of nuclear matter at the high density 
domain by using a description based on the BB model. As in its original form, the 
model is not yet suitable for an extrapolation to higher densities, we considerably 
extent it by considering, at finite temperature, the fundamental baryon octet, lep- 
ton degrees of freedom and trapped neutrinos in a consistent way, i.e. through the 
inclusion of the p-mcson and the additional consideration of the equations for chem- 
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ical equilibrium and charge neutrality. In section 3 we exceed the limit of hadron 
degrees of freedom by including quarks and gluons in order to allow a possible tran- 
sition to quark matter. Within this objective, we focus on a particular process in 
the evolution of compact stars, preceding the supernova explosion: the formation of 
a hot collapsed core or a protoneutron star^ , which can reach temperatures as high 
as few tens of MeV. These hot and dense relativistic objects are formed in a typc-II 
supernova explosion and evolves to a cold neutron star, basically through neutrino 
emission. Here, this very dense and hot core is able to trap neutrinos, imparting 
momentum to the outer layers and then cooling as it reaches a quasi-equilibrium 
state. 

Hitherto, practically all microscopic investigations of stars have been performed 
for dense, cold matter in its ground state. In this work, the evolution of protoneutron 
stars is studied, extending previous investigations to non vanishing temperatures. 
The main novel feature of the present study is the formulation of the dynamical 
equations in the relativistic framework using a mean-field approach. Global static 
properties as masses and radii are then computed as a function of central densities 
and temperatures. 

2. Nuclear Matter 

Nowadays, astrophysical observations of compact stars and terrestrial nuclear 
physics experiments are working together in the determination of the equation of 
state for dense matter. This allows to test the predictions and consistency of differ- 
ent models for the properties of nuclear matter under extreme conditions, through 
the comparison of both astrophysical data and results from high energy heavy ion 
collisions, where nuclear matter approaches similar conditions of density and tem- 
perature as in the interior of dense compact stellar objects as neutron and pro- 
toneutron stars. Conversely, experiments with relativistic heavy ions at current 
facilities are now providing first empirical information on the equation of state at 
high baryon matter densities, particularly relevant to compact stars. Moreover, one 
feature is of particular interest: in addition to the deconfinement regime at very 
high densities and temperatures up to Tc = IbOMeV and beyond, there is another 
phase transition in nuclear matter at much lower excitation energies, analogous to 
the liquid-gas transition in condensed matter physics. This transition is predicted 
to occur with temperatures at Tc > 15 MeV and hence is of particular interest for 
the physics of supernova explosions, and in the formation of protoneutron stars. 
Intermediate energy heavy ion collisions, at energies around 100 MeV per nucleon, 
are currently the most promising tools in searching for the liquid-gas transition in 
laboratory. 

The lagrangian density of our approach, which takes into account the fundamen- 
tal baryon octet, meson sectors, and includes Icpton degrees of freedom, is given 

by 



V- T rl ui , bM ^ c 



4 

4^' 






+ -{d^ad^a - mla^) - -uj^.uj^" + -mluj^uj^ 



This lagrangian density describes a system of eight baryons (B = p, n. A, S^, S'^, 
S+, S~, S") coupled to three mesons (ct, u;, q) and leptons. The scalar and vector 
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coupling constants in the theory, g^ and 17^, and the coefficients b and c for the 
nonhnear a self-couphngs are determined to reproduce, at saturation density po = 
0.153/r7z~^, the binding energy, B = —16.3MeV, and the compression modulus, 
K — 2A0MeV of nuclear matter, and the nucleon effective mass, M* — 732MeV. 
Additionally, the isovector coupling constant gg is determined from the symmetry 
energy coefficient, 04 = 32.5MeV, in nuclear matter. The hypcron/nuclcon coupling 
constant ratios Xi — Qm/gi, with i = cr, w, are constrained through the binding 
energy of the A-hyperon in nuclear matter, from hypernuclear spectroscopy and the 
lower bound of the mass of a neutron star; altogether they determine the ratios of 
the hyperon/nucleon coupling constants, which correlate the relative strengths of 
the a and lo mean-fields to the nucleon and hyperons: 

{B/A)^X^V-XaS, (2) 

and where V = {guI'mi^YpQ and 5* = {M — M*) are the mean-field vector and scalar 
potentials, respectively. 

To determine the nuclear matter equation of state at finite temperature we use 
a set of model parameterizations known as GM3^. With this set of parameters, we 
find for the nucleon effective mass, at saturation density and for the compression 
modulus of nuclear matter, results in the range of current experimental values for 
these quantities as for instance the above cited ones. The effects of nonlinear scalar 
self-couplings are studied by Koepf et. al.^°. Though the ratios Xi ^re known 
only qualitatively, studies have shown however that QHD predictions of bulk static 
properties of neutron stars'* are very sensitive to different choices for Xi- 

Neutrinos, which are not included in the original approach for the GM3 param- 
eters, are an indispensable agency for protoncutron star studies. Neutrinos are thus 
included in our formulation, at chemical equilibrium, via the lepton fraction ratio 

Y, = P^^±J^, (3) 

PB 

where p^, p^ and pb represent, respectively, the electron, neutrino and baryon densi- 
ties. Explicitly, we extract the lepton fraction Yl from the study of the gravitational 
collapse of the core. 

The values for the scalar [a), vector (w) and isovector {g) GM3 coupling con- 
stants, (gi/rrii)'^, with i — a,ui,Q, are shown in table 1 together with the hy- 
peron/nucleon coupling constant ratios, and the lepton fraction, Y^. 



Table 1: Parameters adopted in the model. 



{9a/m„Y {g^/m^Y {gJnigY b (xlOO) c(xlOO) 

fm^ fw? fra^ 

9.927 4.820 4.791 0.8659 -0.2421 

Xa Xlo Xq Yl 

0.6 0.568 1.0 0.4 



Applying standard technics from field theory and the mean-field approach, we 
obtain the equation of state for nuclear matter as a parametric equation p = p{e) 
which explicitly relates the energy density 
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2Jb + 1 



E 

B 

E 



(27r)3 
2 

WrJo 



^k^ + Mf[n{fiB, T) + n(MB, r)]ci3fc 



v/fc2 + mA2[n(^A, T) + n(AiA, r)]d3fc ^ 



(4) 



and pressure 
P = 



2 2 

-m„(7 



1 



-E 



2 ' 
2Jb + 1 

(27r)3 
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2 2 
fc2 



2™e£'03 



fc2 + M^2 



6M 3 c 4 
3 4 

[n(fj,B,T) + n{pB,T)]d^k 



3 ^ (2^)3 ^0 



\/fc^ + JTiA^ 



[n{iix,T) + n{^ix,T)]d^k. 



(5) 



In the expressions above we may identify, in the first row of both, the contributions 
of the scalar, vector and isovector mesons, followed by the contributions of the 
Fermi gas of baryons and leptons; these contributions appear as sums over the 
particle families, B = n,p, I],S,A and A = e,^^v, and where Jb is the isospin 
baryon number. The quantities n(/x, T) and n(/i, T) are the Fermi-Dirac distribution 
functions for particle and anti-particles, at finite temperature: 



n{p.B,T) = [e' 



/3(^k2+m-2-^b + ^Pb + ^P3) 



ir\ 



n{^iB,T)^[e - e +1] , 



(6) 



with (3 = 1/kT and the chemical potentials ^b and ^\ for baryons and leptons, 
respectively (see Figs.l and 2 where the Fermi distribution function and the density 
of states are plotted as functions of chemical potential and the squared momentum 
e, respectively). 
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Figure 1: Fermi distribution rife as a function of the baryon chemical poten- 
tial for T=OMeV (solid line), T==10MeV (dotted line), T=50MeV (dashed hne), 
T==100MeV (long dashed hne) and T = 200MeV (dot-dashed line). 
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Figure 2: Density of states as a function of baryon squared momentum for T=OMeV 
(solid line), T=50MeV and T^lOOMeV 



In the numerical calculations with the BB model, we have considered hyperonic 
matter without (Fig. 3) and with trapped neutrinos (Fig. 4), and at a higher tem- 
perature (T=50MeV in Fig. 5) in order to evaluate the effect of these quantities on 
protoneutron star properties. 




Pb (fm " ) 

Figure 3: Distribution of particles as a function of baryon density for T=10MeV, 
without neutrino degrees of freedom. 

A liquid-gas phase transition is predicted in various relativistic nuclear mod- 
els such as QHD, BB and others, when symmetric nuclear matter or pure neutron 
matter are described. However, following within the same models to describe asym- 
metric nuclear matter, via the inclusion of the fp-mcson, this phase transition disap- 
pears. Instead, the energy used to generate this liquid-gas phase transition, leads 
now to hyperon production, resulting in a phase transition of asymmetric nuclear 
matter into hyperonic matter. In such a phase transition, the system has to provide 
an extra energy to convert nucleons into hyperons which makes the corresponding 
EOS stiffer. The provision of this extra energy is not easily seen in this case in 
the EOS (Fig. 6), as it is in the former Van der Waals liquid-gas phase transition. 
It is however clearly visible in the curve for the specific heat, where each hyperon 
generation formed introduces a discontinuity on its structure (Fig. 7). This phase 
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Pa (fm ) 

Figure 4: As figure 3, including however neutrinos in the calculations. 




Figure 5: As figure 4, however for T=50MeV. 



transition is similar to the corresponding one with the formation of Cooper pairs in 
the supercondutivity phenomena in electron systems, where discontinuities are also 
found instead of singularities, characterizing this type of behaviour as a continuous 
phase transition. 

3. Quark Matter 

We have developed so far, nuclear matter studies disregarding the presence of 
quarks and gluons degrees of freedom. However, the extreme conditions for the 
density in the interior of a protoneutron star may be favorable for the occurrence of 
a hadron-quark phase transition, which may form a core of a quark-gluon plasma 
(QGP) in the inner regions of the star. 

We make use of the MIT bag model^^ in order to describe the quark matter 
equation of state and the hadron-quark phase transition. The MIT lagrangian is 
given as 



Cmit = [k^l'^df.^b - (d^^^h^i;) - B\Qv(x) - \^^^s , 



(7) 



(for details on the meaning of the symbols, see for instance^^). Following again stan 
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Log,,, e (g/cm) 



Figure 6: Nuclear matter equation of state for T=OMeV (solid line), T=50MeV 
(dotted line) and T=100MeV (dashed line). 



P (fm ) 



Figure 7: Specific heat as a function of baryon density for T=10MeV. 



dard technical procedures of quantum field theory, this lagrangian density generates 
the following equation of state for quark matter {p ~ p{£))' 



f 2" ^° 



(8) 



p = -B + Y,^^ J fc'rffcK(M)+n/c(A*)] 



(9) 



with the degeneracy factor represented by 7/ = 2spin x Scoiour and the bag constant 
denoted by B. The baryon density is given by 



..-E^^f'^^M 



nk{^J.) -rikip)]. 



(10) 
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The hadron-quark phase transition is determined via the Gibbs criteria^^, which 
require that, at constant 

temperature, pressure and chemical potential of both hadron (H) and quark (Q) 
phases are related by 



Pq = Ph, 

for a conserved overall baryon number, i.e. 

Nq/3 + Nh ~ constant , 



(11) 



(12) 



where Nq is the number of quarks and Nh is the number of baryons in the system. 
The resulting equation of state includes the presence of hadrons, with a phase 
of deconfined quark matter represented in Fig. 8, and where the plateau shown in 
the figure reflect the Gibbs condition of the phase transition. 



-9.5 
Log,,, e (MeV/fm ') 



Figure 8: Equation of state with the hadron-quark phase transition for a bag 
constant B = 100 M eV / fm^ (solid hne), B = l3l.2MeV/fm^ (dotted line) and 
B = IbOMeV/fm^ (dashed line). 



4. Stellar Structure 

The protoneutron star configuration may be determined by combining the nu- 
clear matter EOS with the TOV equations from general relativity: 



Airr dp(r) ~ — 



M{r)dm(r) 



PJr) 
e{r) 



M{r) 



2M(r) 



(13) 



and 



M{r) 



47r / eir'y^dr'. 

10 



(14) 



The above equations describe the structure of a static, spherical and isotropic 
star wit the pressure p{r) and the energy density e(r) reflecting the underlying nu- 
clear model; M(r) denotes the mass inside a sphere of radius r. The TOV equations 
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involve various constraints and boundary conditions: they must be evaluated for 
the initial condition e(0) = ec (with ec denoting the central density) and M(0) = 
at r = 0; the radius R of the star is determined under the condition that on its 
surface the pressure vanishes, p{r)\r=R = 0. 

The condition for chemical equilibrium for protoncutron stars are: 



A^i = hy-n -qi{l^£ - l^i^e) 



(15) 



where fXi and iii stand for the baryon and Icpton chemical potentials, respectively; 
bi is the baryon number and the baryon and Icpton electrical charges are represented 
by gi- ^ ^ 

The corresponding equations for baryon number and electric charge conservation 
are: 



Pb. 



l3 
V^ '^F,B 



and 



Y.1- 



'aryomc 



^F.B 



37r^ 



^3.2 



37r2 



(16) 



(17) 



In the figures below, the bulk properties of the protoncutron star as a function 
of the temperature are shown. Figs. 9 and 10 show respectively the maximum mass 
and its radius for hyperon matter including neutrino trapping (solid curve) and 
the inverse situation, where neutrinos escape freely (dotted curve) These curves 
indicate maximum masses respectively equal to Mum = 1.547Mq for the former 
case and Mum = 1.557M0 for the latter. Inspection of these figures indicates that 
for temperatures higher than 70 MeV equilibrium is obtained for lower masses and 
radii. If T < 70 MeV the limiting mass increases with temperature. 



50.0 

T (MeV) 



Figure 9: Protoncutron star maximum mass as a function of temperature for Yi 
0.4 (solid line) and Yl = (dotted hne). 



The analysis of the maximum mass vs. temperature relationship is not trivial 
because of the complex structure of the TOV equations and those describing the 
thermodynamical state of matter. In one hand, to the increase of temperature it 
corresponds an increase of the mass of the protoncutron stars, a very well known 
result for the model used in this calculation*; however, other models predict exactly 
the opposite behaviour. The novel feature we have found without any doubt is an 
absolute limiting value for the protoncutron star mass, followed by a subsequent 
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50.0 

■ (MeV) 



Figure 10: Radii of protoneutron stars with niaxinium-niass as a function of temper- 
ature for trapped (solid line) and for free (dotted line) neutrino matter, respectively. 



decrease of this maximum value with the increasing of temperature. The existence 
of this maximum mass for protoneutron stars, reached at T ^ lOMeV is a new result 
obtained in this study and we mention initially three possible reasons to explain 
it: the very sensible balance of pressure, energy density and gravity in the TOV 
equations; the modification of the lower limit on the energy density determined 
by thermal contributions, and the complex behaviour of the scalar field at finite 
temperature. Here we briefly comment on the last two aspects and in the following 
we present some analytical estimates for this behaviour. 

As temperature increases, due to thermal contributions, an increasingly minimal 
value for the energy density of nuclear matter is progressively reached. Our results 
indicate that, when taking into account into the calculations this characteristic 
behaviour of the thermal energy density, it implies sharper crusts for T > 0, as 
compared to the corresponding results at zero temperature. As this effect starts 
to be significant around T — 50MeV, it might be a signal that the maximum 
protoneutron star mass starts to decrease. 

A more microscopic and most interesting explanation for the results previously 
stressed, with striking consequences for nuclear matter at high densities, can be 
found in the opposite behaviour of the meson scalar field at the high and low density 
regimes. At low baryon densities, the scalar field does not change significantly as a 
function of temperature until T reaches several tens of MeV; then, from this value 
on, it starts to increase. On the contrary, for high baryon densities, the scalar meson 
field decreases as the temperature increases already beyond the first few MeV. As 
the temperature increases, this behaviour becomes even more significant and thus 
might be basically responsible for the protoneutron star maximum mass behaviour 
exhibited here. 

The temperature dependence of the maximum star mass is supported by ana- 
lytical estimates of the T-dependence of the compression modulus and the equation 
of state of nuclear matter. 

Concerning the compression modulus of nuclear matter, 



Kipo,T) 



9p§ 



dp 
dp 



'P=PO 
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9p§ 



dp dp f dp 

'dp '''=''" ^ 9r \df 



'P=PO 



we find for its leading dependence in T 

K{po,T) ^ K{po,T ^ 0) (1 + aT^) , 



(18) 



(19) 



with a — 0.0181. This expression indicates that as temperature increases, the 
compression modulus of the protoneutron star matter also increases, giving a stiffer 
equation of state, for a fixed density, in comparison to the case T = 0, with the 
immediate consequence that the stiffer the equation of state, the more massive 
becomes the limiting star mass of a family of compact stars. 

A second analytical argument involves the threshold equation for a given species 
of hadrons 



qsp-n 



qe,BPe > gLuB^O + 9gBQidzhB + [Mb - 9aB(^) ■ 



(20) 



where the isospin term ggBQmluB in the eigenvalue determines whether a species is 
isospin-favored or unfavored. By using for the chemical potential^^ 



with 



i^B{T 



Pb-^Pb{T = 0) 


1 / ^kT V 




3 \2pb{T = 0)) _ 


7 


3) = giuB^O + ggBQmhB + V fc| + (Mb 


- gaB(7f 



(21) 



(22) 



one finds (see also figures (3)-(5)), that at densities relevant for neutron stars the 
relative density distribution of the hyperons, which contribute to the softening of 
the equation of state and thus to a reduction of the limiting star mass, is shifted 
to higher values with increasing temperature. As a consequence, the limiting mass 
of the protoneutron star will increase in the comparison with the corresponding 
value at T = 0. Similar arguments were recently used by Prakash^ in his studies on 
compact stars a finite T. 

Apart from these considerations, a more complete analysis should involve both 
TOV equations which allow, for an appropriate discussion, only numerical solutions. 
On the other hand, still some insights may be obtained from the T-depcndence of 
the energy density e(K,T). Taking the Fourier transform at finite temperature, 



e(r,r) = 



1 



(27r)3/2 
1 



e{K,T)e'^-''d'K 



(2^)3/2 7o 'Jo 'Jo 



4:71 



(2^)3/2 7o 



dn K jo{Kr)e{K,T) . 



(23) 



and combining this expression with equation (14), we find the following equation 
for the protoneutron star mass 



M ■ 



AttR 



dKe{K,T^O)ji{KR) + 



AttR 

(2^2 7o 



dKes{K,T)j,{KR). (24) 
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Above, R represents the protoneutron star radius and k is the Fermi momentum; 
furthermore, ji{KR) with i = 0, 1 represent Bessel functions of i-kind while 65 (k, T) 
denotes the temperature corrections described in the appendix 

es(At,r) = 5]-|(fcBT)2[2fcf,,i?f,, + -^] (25) 

with kp^B — K (see equation (A. 5)). Evidently, one identifies the first term in 
expression (24) as the protoneutron star mass at T=0 and the second term as 
finite temperature corrections; the Bessel function in expression (24) modulates the 
kernel of the integral. As a crucial point, from the spectrum of the energy density, 
low momenta k (with K.R below the first zero of the Bessel function) dominate the 
integration over k and add up coherently, resulting in an increase of the protoneutron 
star mass as a function of T. 

We remark that our conjecture is supported both by a quite similar result in 
the non-relativistic limit obtained by Chandrasekhar^^ in the description of the so 
called Chandrasekhar mass (see also Baron and Cooperstein^^) and from the studies 
recently performed by Woosley and Weaver on the birth of neutron stars^^, who 
find for the temperature correction in the protoneutron star mass 



M«Mt=o(1 + ^2 ) (26) 

with ep being the Fermi energy. 

The inclusion of quarks has a significant impact on the maximum mass of a 
neutron star: hybrid stars may develop a QGP core, which may extend several 
kilometers, lowering its maximum mass due to the stiffening of the EOS. For this 
reason, to develop a treatment compatible with the increasing of this value, we have 
adopted here different parameterizations for the nuclear matter equation of state. 
At even higher central densities one can find other classes of stars, for instance 
strange and quark stars, which are basically formed by deconfined quark matter 
with a thin (1 km) crust of nuclear matter. 

However, there is a very important condition which shall be taken into account 
when phase transitions are considered in neutron stars. The Seidov criterium^'' 
establishes that whenever a density or/and an energy density jump exists inside the 
star, consequence of a phase transition of any type, the resulting configuration is 
stable only if satisfies the condition 

ei 2 V ei/ 

where p is the pressure at the interface between both phases, while e2 and ei are 
respectively the energy densities for the higher and lower density phases. It should 
be emphasized that this criterium is applicable for fully relativistic star models. 

The models described in this work generate hybrid proto-neutron stars which 
satisfy that condition at zero temperature. However, objects hotter than 20 MeV 
are not able to satisfy the Seidov criterium and become unstable as soon as they 
underwent a phase transition. 

5. Summary and Conclusions 
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In the framework of the Sommcrfcld approximation, we have studied the struc- 
ture of protoneutron stars through an effective many-body field theory at finite 
temperature. To the best of our knowledge, the approach used in this work, includ- 
ing the Seidov criterium, was followed for the first time in the present investigation 
of this topic. The Sommerfeld approximation permits a drastic simplification of 
the computational input and thus provides a rather clear insight in the role of 
temperature on bulk static properties of these fascinating stellar objects. 

The calculations were performed by using an extended version of the non-linear 
BB model, including the fundamental baryon octet, a, uj and g mesons, leptons 
degrees of freedom, trapped neutrinos, — introduced into the formalism at fixed 
lepton fractions, — the equations for chemical equilibrium and charge neutrality. 
Integrating the TOV equations we have obtained the maximum mass and the cor- 
responding radius of the protoneutron star as functions of the temperature. 

The results indicate that to the increasing of temperature, it corresponds the 
increasing of the maximum mass and radius of the star and, moreover, that the 
EOS of nuclear matter becomes stiffer. The causal limit is exactly recovered in the 
EOS, even in the Sommerfeld approximation, where its recovering becomes more 
transparent. However, at very high temperatures, the EOS acquires, compared to 
the T = results, a lower limit on the energy density which prevents the increasing 
of the maximum mass of the protoneutron star; beyond a critical temperature, 
which we have found equal to T = TQMeV , the maximum mass of protoneutron 
stars decreases. 

We have analyzed the protoneutron star structure with and without neutrino 
trapping and the corresponding maximum masses obtained are 1.557 Mq and 1.547 
Mq, respectively, for the set of parameters shown in table 1. We have identified 
this value as the maximum protoneutron star mass which can be formed right after 
the supernova event. With the same parameters in table 1, except for the choice 
of the hyperon-nuclcon coupling ratios^ Xa = Xi^ — Xg ^ ■\/2/3, we have obtained 
respectively 1.889 Mq and 1.868 Mq. These last values are able to explain the mass 
of the neutron star in the binary system Vela X-1^ and this example illustrates how 
these mass determinations may guide us in the choice of adequate nuclear interaction 
coupling constants. 

In figure 11 we exhibit the behaviour of the protoneutron star mass as the 
temperature increases (in order to group three different stars in the same plot, we 
normalize the protoneutron star masses by their values at zero temperature) . Each 
of the three sets of stars are chosen to have the same total baryon number. 

We also have seen that a hybrid protoneutron star, which shall obey the Seidov 
criterium, can only be formed at low temperatures of the order of T=20MeV. Pro- 
toneutron stars hotter than that value shall first cool down as hyperon stars and 
may form, subsequently, a quark core in its inner region. 

The density dependence of the scalar meson field is of particular interest for dense 
hadronic systems in general, as it touches upon the basic question how to organize 
effective field theories. Here two scenarios are obvious, which lead to a drastically 
different dependence of the nonlinear scalar (and vector) meson field with density 
and temperature. In the naive dimensional analysis, the scale parameter for natural 
nonlinear self-coupling of the a and io fields is the effective mass scale 

meff -- ^/JJv - (250 - 300MeF) , 

at saturation density, with /^ ^ ^iMeV as the weak 7r-decay constant and A ^ 
nucleon mass as the scale for chiral symmetry breaking. From arguments based 
on Brown-Rho scaling and the onset of deconfinement with the appearance of a 
quark-gluon plasma beyond some critical density pc ^ 4:pQ, both /jr and A and thus 
rueff decreases with increasing p yielding 

meff ^Oforp>pc. 
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T (MeV) 

Figure 11: Protoneutron star mass normalized by its mass at T=0, as function of 
temperature for Mt=o = 1-SMq (solid line), Mt=o = 1-6Mq (dotted line) and 
Mt=o = 1.35M0 (dashed line). 



Thus organizing effective field theory with the scale mass -^//jrA would completely 
change the role of high order meson self-couplings — in the limit vrieff -^ the 
corresponding lagrangian would develop an essential singularity — and at the same 
time, used as a scale parameter for the temperature dependence 



A 



fTT-eff 



would invalidate the Sommerfeld approximation beyond a critical density. 

However, an alternative scale for fermionic matter is set by the Fermi-momentum 
kp itself. At saturation density the two scales coincide: at p = po 

270 MeV . 

Opposite to rrieff, however, kp is stable, as it increases typically as kp ^ p^'^ with 
increasing density. Such an organizational scheme of effective field theory would 
preserve both the dominance of the cubic and quartic a (and lo) self-couplings and 
with the expansion parameter 



A 



kT 



kT 



kp (270Mey)(p/po)i/^ ' 



would even improve the convergence of the Sommerfeld approximation with 
increasing density. 

Our investigations are just a first attempt towards the comprehension of the 
role of temperature on physical properties of neutron stars. Our main findings are 
substantiated by calculations based on improved microscopic models, which include 
as an example, in the spirit of the Brown- Rho scaling, a density dependence of the 
coupling constants in the model. Investigations along these lines are in progress. 

Appendix A 

The finite temperature studies involves many technical difficulties through the 
Fermi integrals; in particular, the computational effort increases dramatically. How- 
ever, as from the integration of the full energy density only a narrow range con- 
tributes significantly to observables, approximate schemes should provide at least 
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semi-quantitative results. The Soinmerfeld approximation is a very successful ap- 
proach first applied in condensed matter physics and recently applied to the white 
dwarf problem^^. We briefly describe, in this section, the Sommerfeld approxima- 
tion and discuss its validity when applied to nuclear matter problems. The main 
conception of the approximation is the following: if a given function H{e) of the en- 
ergy density e does not vary rapidly in the energy range of the order of kT around /i, 
the temperature dependence of the integral /_ H{e)n{e)de should be given quite 
accurately by replacing H{e) by the first few terms in its Taylor expansion around 



°° in 
n=0 



(e-/z)" 



e=M 



The result is a series of the form 

/oo r-fj. °° j2n— 1 

H{e)nie)de = / H{e)de + Y^{nTf-a^-—-^H{e) 



(A.l) 



(A.2) 



The coefficients a„ are dimensionless constants of order 0(1), and are determined 
according to 

' dx. (A.3) 



, (2n) V dxe^ + 1 

The functions H{e) of our approach have major variations on energy scales of 
the order of /x and thus in general ((i/(ie)"i/(e)|(,=^ is of the order of H{ii)/fi". 
Consequently, successive terms in the Sommerfeld approximation are smaller by 
0((kT)^) ~ 0(10"^). Consequently, in actual calculations only the first few terms 
must be retained. The final form of the integrals, up to the third order, is then 
given as 

H{e)f{e)de = /'' H {e)de + ^ {nTf H' (fi) 

_j J — oo ^ 

+ !^(^TrH"'ip)+Oi — )\ (A.4) 

360 yU 

where we recognize the zero temperature contribution corrected by terms which 
include explicitly the temperature. 

We describe, in the following, the EOS of nuclear matter at finite temperature, 
extending the T = energy density eq and pressure po, by the leading T-dependent 
coefficients from the Sommerfeld approximation: 3 

« = £0+ V ^(A:BT)2[2fc^,,i?^,, + -^] (A.5) 



3 / n \ 5 



^''<^«-''l-|f-'(tf) "(If 



Po 



2i(,,rnf£^-fi| (A^e, 



i=B, 



7 2n rT,\A,rr,A^F,i on I ^F,i\ , n I -^F,i\ , ^^(Epj 



72d''"^=^''P^^-^'liJ?J +n^J +'\«« 
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In previous appreciations, the Sommerfeld approximation has been used to de- 
scribe systems in the high density and low temperature ranges, tt^ >> 1, a con- 



ksT ^ 

dition which is very weU satisfied in protoneutron stars, with typicaUy jr^ 



[10 — 10^]; with this in mind, we feel confident about the adequacy of the approach 
for the investigation here exhibited. However, to explore quantitatively the limiting 
range of validity of the Sommerfeld approximation, we have performed calculations 
for less complex microscopic models and compared the approximate results with 
numerically rigorous calculations (sec figure (A.l). 




Logj,, e (MeV/lm ) 

Figure A.l: EOS for the Walccka model for the temperatures T = 50 MeV and 10 
MeV (upper and lower curves, respectively). Solid lines: exact results; dashed and 
dotted lines: first and second order corrections of the Sommerfeld approximation. 

We have used, in particular, the Sommerfeld approximation to describe the 
equation of state of the original QHD model for pure neutron matter. Then, the 
equation of state was combined with the TOV equations, giving the results shown 
in table A.l for the exact calculations and the approximated ones using the Som- 
merfeld based calculations. We have tested our results at the critical temperature 
for the liquid-gas phase transition in nuclear matter (T=9.2MeV) and for higher 
values (T=50MeV), as typically found in protoneutron stars. For comparison, the 
maximum mass predicted in the Walecka model at T = is 2.57 Mq. 



Temperature Exact Sommerfeld 



9.2MeV 
50MeV 



2.606 
3.215 



2.605 
3.144 



Table A.l: Comparative results of the maximum protoneutron star mass. 



The (nearly) quantitative results of the two approaches confirms the adequacy 
of the utilization of the Sommerfeld approximation in the description of physical 
properties of nuclear matter at finite temperature. 
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